Controlling Luttinger Liquid Physics in Spin Ladders under a Magnetic Field 
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We present a ^^N nuclear magnetic resonance study of a single crystal of CuBr4(C5Hi2N)2 (BPCB) 
consisting of weakly coupled spin-1/2 Heisenberg antiferromagnetic ladders. Treating ladders in the 
gapless phase as Luttinger liquids, we are able to fully account for (i) the magnetic field dependence 
of the nuclear spin-lattice relaxation rate Tj~^ at 250 mK and for (ii) the phase transition to a 3D 
ordered phase occuring below 110 mK due to weak interladder exchange coupling. BPCB is thus an 
excellent model system where the possibility to control Luttinger liquid parameters in a continuous 
manner is demonstrated and Luttinger liquid model tested in detail over the whole fermion band. 

PACS numbers: 75.10.Jm, 75.40.Cx, 76.60.-k 



Interaction between quantum particles plays a crucial 
role in one dimension (ID), where its interplay with quan- 
tum fluctuations leads to a state described as a Luttinger 
liquid (LL) [l|. Low-energy physics of the LL is fully 
characterized by two interaction dependent LL parame- 
ters: the velocity of excitations, u, and the dimensionless 
exponent K. Correlation functions decay as power laws, 
with exponents which are simple functions of K. The 
LL model has been shown to apply to a growing number 
of ID systems, such as organic conductors 0, quantum 
wires [sl, carbon nanotubes [3], edge states of quantum 
Hall effect [5] , ultra cold atoms [6| , and antiferromagnetic 
(AFM) spin chain 7] or spin ladder systems Several 
characteristic features of the LL model have been ob- 
served in these systems, such as the power law behavior 
of some correlation or spectral functions. However, since 
the details of the interaction are rarely known, only a the- 
oretical estimate of the power law exponents is usually 
possible. A precise quantitative check of the LL model is 
thus still missing. 

The obstacle can be overcome in spin ladder systems. 
Namely, spin-1/2 AFM ladder in an external magnetic 
field H maps essentially onto a ID system of interact- 
ing spinless fermions [1, 9, [13, 11|, where H acts as a 
chemical potential. The interaction term in the fermion 
picture is uniquely determined by the exchange coupling 
constants, which can be experimentally extracted [J± on 
the rungs and J|| on the legs of the ladder; see Fig. 1(a)], 
and by H, which controls the filling of the fermion band. 
By increasing H the spin gap between the singlet ground 
state and the lowest triplet {Sz = — 1) excited states of 
the spin ladder decreases. It closes at Hd, where the 
ladder enters the gapless phase corresponding to the par- 
tially filled fermion band. At Hc2 the ladder gets fully 
polarized and the gap reopens. Once J± and J|| are 



known the LL parameters can be obtained numerically 

for arbitrary H. The associated LL prediction can then 

be checked quantitatively over the whole fermion band, 

which extends between the critical fields Hd and Hc2 ■ 
Suitable spin ladder systems are rare, either because 

of unattainable critical fields or because of the presence 
of anisotropic interactions, such as Dzyaloshinski-Moriya 
(DM) interaction [13], as is the case in extensively stud- 
ied Cu2(C5Hi2N2)2Cl4 Recently, CuBr4(C5Hi2N)2 
(BPCB) [ij] was identified as a good candidate. Namely, 
low-temperature magnetization data were well described 
by the XX Z chain mod el [isl in the strong-coupling limit 
{Jj_ ;2> J||) of a ladder loT IsL 13] • In addition, the re- 
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FIG. 1: (color online), (a) A ladder formed by the supposed 
exchange interactions, JijSi ■ Sj, between S — 1/2 spins of 
Cu^^ ions in the crystal structure of BPCB. Solid and dashed 
thick blue lines stand for J± and Jy , respectively. The 10 
protons attached to the C atoms are not shown, (b) ^^N 
NMR spectra at 120 and 40 mK recorded at = 9.0 T. 
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suits of thermal expansion and magnetostriction experi- 
ments were explained within the free fermion model and 
in more detail with quantum Monte Carlo calculations for 
spin ladders 17, 1^. In this Letter we present a detailed 
^^N nuclear magnetic resonance (NMR) study of BPCB 
revealing for the first time the occurence of the field in- 
duced 3D magnetic order below 110 mK. We show that 
the LL model completely accounts for the experimental 
behavior of BPCB in the gapless phase. The determined 
phase diagram and field variation of the order parameter 
are perfectly described in the framework of weakly cou- 
pled LLs. The importance of these results is twofold: (i) 
they show that the whole physics of coupled spin ladders 
can be captured in a single theory based on LLs, and (ii) 
they provide the first quantitative check of the LL model. 

In ^"^N NMR experiments we used a single crystal of 
BPCB with dimensions 3.5x2x4 mm^. As shown in 
Fig. 1(a), the pairs of spin-1/2 Cu^+ ions (spin dimers) in 
the crystal structure of BPCB are stacked along the crys- 
tallographic a axis to form ladders [13] ■ Each unit cell 
contains two rungs of two different ladders (denoted by I 
and II), which are crystallographically equivalent, but be- 
come physically inequivalent when the external magnetic 
field H is applied in an arbitrary direction. There are 
only two crystallographically inequivalent nitrogen (N) 
sites per ladder in a unit cell: N(l) is located close to the 
rung and N(ll) close to the leg of the ladder [Fig. 1(a)]. 
Since "'^^N has spin 1=1 and thus a quadrupole moment, 
each site gives rise to a pair of ^''N NMR lines (doublet), 
split by the quadrupole coupling with the local environ- 
ment. Therefore, -'^^N NMR spectrum [upper spectrum in 
Fig. 1(b)] consists of four doublets, for two sites in each 
of the two ladders. We oriented the sample so that the 
external magnetic field lay in the a*b plane (a* _L b,c), 
at an angle of 9° to the b axis. In this orientation, the 
N(l) doublets are well resolved, while the N(ll) doublets 
overlap [Fig. 1(b)]. All the measurements presented here 
were performed on the N(l)i lines. 

The magnetization of the Cu^+ ions is detected via 
the hyperfine coupling to the ^"^N nuclei. It is easily re- 
constructed due to the inversion center in the middle of 
each rung [3] ■ Any uniform magnetization of Cu^"*" ions 
shifts the average position of the NMR quadrupole dou- 
blet with respect to the Larmor frequency. Any staggered 
(i.e., AFM) magnetization of Cu^+ ions, however, breaks 
the inversion symmetry between the pairs of equivalent 
N sites. This doubles the number of inequivalent N sites 
and each NMR line splits symmetrically in two lines. The 
associated hyperfine shift and splitting are proportional 
to the uniform and staggered magnetization of Cu'^^ ions, 
respectively. Fig. 1(b) illustrates the emergence of the 
staggered magnetization in BPCB aX H = 9.0 T on cool- 
ing from 120 to 40 mK. Its origin is discussed later. 

Next we check whether our longitudinal (parallel to H) 
uniform magnetization mz{H) data are well described in 
the XX Z chain model. Fig. 2(a) shows mz{H) (per spin) 



measured at 40 mK via the hyperfine shift of N(l)i NMR 
lines, from which we extract the values of the critical 
fields: ffci = 6.703±0.008T andiJc2 = 13.888±0.006 T. 
Instead of the square-root singularities expected for a 
perfect ID system, very close to the critical fields we ob- 
serve a linear mz{H) dependence [Fig. 2(a) inset], a sig- 
nature of weak inter-ladder (3D) exchange coupling [lO|. 
As this 3D regime is very narrow, by fitting mz (H) data 
with a phenomenological function having square-root sin- 
gularities at ?7iz = and 0.5 [Fig. 2(a) inset], we can 
estimate the values of the corrected, "ID" critical fields 
pertaining to an isolated ladder: H^^ = 6.763 T and 
H^^ = 13.828 T. With g = 2.176 corresponding to our 
sample orientation 14] this yields [l^ J± = 12.9 K and 
J|| = 3.6 K, hence JjJ J\\ = 3.6, in agreement with recent 
determination IT], ll8| . With these values we perform the 
density matrix renormalization group (DMRG) calcula- 
tion of niz (H) at T = for a single ladder and find an ex- 
cellent agreement with our experimental data. As shown 
in Figs. 2(a,b), the slight asymmetry of the curve about 
its mid-point at mz = 0.25 is nicely reproduced since 
we take into account also the upper two triplet states of 
each dimer {Sz = 0, 1). These are neglected in the strong- 
coupling treatment {J± ^ J|| ) leading to the XX Z chain 
model [lO| . The corresponding curve is symmetric about 
its mid-point, in worse agreement with our experimental 
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FIG. 2: (color online), (a) Magnetic field dependence of 
the longitudinal uniform magnetization per Cu^'*' ion at 
40 mK (o). Inset shows linear mz(H) dependence very close 
to Hc2 (solid line) and the determination of (dashed 
line) . The data are compared to the result of the DMRG cal- 
culation for J±/J\\ = 3.6 (solid line) and to the prediction of 
the XXZ chain model (dashed line), both at T = 0. (b) Dif- 
ference between each prediction and the experimental data, 
(c) Variation of the LL parameters K{H) and u{H) (in kelvin 
units) over the fermion band as calculated for J±/J\\ = 3.6 
(solid line) and for the XXZ chain model (dashed line). 
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data. However, mz{H) as a thermodynamic quantity is 
not very sensitive to the model. Once the couphng ra- 
tio J±/J\\ is fixed, the variation of the LL parameters 
with H is completely determined and the LL theory is 
left without any adjustable parameter. For an isolated 
ladder with J±/J\\ = 3.6 we numerically calculate K[H) 
and u{H) (in kelvin units), combining DMRG method 
with bosonization as in Ref. [lH. The result is displayed 
in Fig. 2(c) together with the corresponding result in the 
XX Z chain model (from Ref. 0). Close to Hlf [Hlf) 
the LL exponent approaches the value K = 1 oi nonin- 
teracting fermion system indicating a nearly empty (full) 
fermion band. 

The LL behavior is tested via the dynamical spin-spin 
correlation functions, which are experimentally accessible 
through NMR observables. We focus on three such ob- 
servables, starting with the nuclear spin-lattice relaxation 
rate T^^ . At low temperatures T^^ probes exclusively 
Cu^"*" spin dynamics, namely its low-energy (w — > 0) part 
[iol corresponding to the long time behavior of the local 
spin-spin correlation functions. For the above determined 
range of K{H), by far the biggest contribution to Tj~^ at 
low temperature comes from the transverse (perpendicu- 
lar to H) staggered correlation , and we find it to be 
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where Aq is the amplitude of the correlation function, 
j/{2tt) — 3.076 MHz/T is ^^N nuclear gyromagnetic 
ratio, A±_ the transverse hyperfine coupling constant, 
and B{x,y) = T{x)r{y)/T{x + y). Fig. 3(a) shows 
T^^{H) dependence measured in the gapless phase of 
BPCB at 250 mK, well above the 3D ordering tempera- 
ture. Its concave shape reflects the increased ID fermion 
density of states close to the critical fields. Using the 
previously calculated K{H), u{H) [from Fig. 2(c)] and 
Aq{H) (calculated along with K and u) for J±/J\\ — 3.6, 
T^^{H) from Eq. ([T|) is compared to the data by ad- 
justing a single scaling factor A"^ [Fig. 3(a)]. Very good 
agreement over the whole field range provides a remark- 
able confirmation of the LL model. Moreover, the utilized 
value A± = 570 G agrees with that obtained from direct 
-'^■*N NMR determination [2l[. In contrast, the curve ob- 
tained in the XX Z chain model (with the same A±) fails 
to reproduce the biased shape of rj"^(iJ) [Fig. 3(a)]. This 
demonstrates the sensitivity of the observable T^^{H) to 
the applied set of LL parameters K{H), u{H) and, con- 
sequently, to the coupling ratio J±/J||. 

The (dominant) transverse staggered spin-spin cor- 
relation function diverges with decreasing temperature, 
which leads at low enough temperature to a 3D ordering 
due to weak interladder exchange coupling [l3] . The as- 
sociated order parameter is a transverse staggered mag- 
netization nix (per spin) , measured via the hyperfine 
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FIG. 3: (color online). Magnetic field dependence of (a) 
Tj"^ at 250 mK (o), in the ID LL phase, (b) the temperature 
of the transition between the 3D ordered phase and the ID 
LL phase (o), (c) the transverse staggered magnetization rux 
per Cu'^"'' ion at 40 mK (o) and its extrapolation to zero 
temperature (•), all measured on N(l)i site as described in 
the text. Each data set is compared to the corresponding 
prediction of the LL model based on K{H) and u[H) [from 
Fig. 2(c)] for J_l/ J\\ = 3.6 (solid red line) and to the prediction 
of the simplified XXZ chain model (dashed green line). The 
only adjustable parameter in each case is an overall scaling 
factor. Dotted blue line in (c) is the result of the DMRG 
calculation. 



splitting of the N(l)i NMR lines. To map the boundary 
between the 3D ordered phase and the ID LL phase, we 
determine the temperatures Tc{H), at which this split- 
ting vanishes [Fig. 3(b)]. Since the whole phase boundary 
lies below 110 mK (^ J|| = 3.6 K), we expect that the in- 
terladder coupling (i.e., the coupling between the LLs) is 
small and treat it in the mean-field approximation lead- 
ing to 0: 
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Here J' is the exchange coupling between the Cu^+ ions 
in neighboring ladders and z the coordination number. 
Although exchange paths between the ladders have not 
been identified yet, the structure of BPCB suggests that 
z^A. Using K{H), u{H) and Ag(iJ) for J_l/J|| = 3.6, 
we compare Tc{H) from Eq. ((2|) to the experimental data 
by adjusting a single scaling parameter J' [Fig. 3(b)]. For 
J' = 20 mK [2^ we obtain an excellent agreement, with 
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faithfully reproduced biased shape. The prediction of the 
XX Z chain model (with the same J') is considerably dif- 
ferent in shape, showing that Tc{H) is another observable 
very sensitive to the applied set of LL parameters. 

The field dependence of the order parameter mx{H) is 
measured via the hyperfine splitting of the N(l)i NMR 
lines at 40 mK [open circles in Fig. 3(c)]. The data are 
compared to the T = prediction, which we obtain in 
the mean-field treatment based on weakly coupled LLs: 



= F{K) 

where, using the results from Ref. 24| . 
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Taking again K{H), u{H), A^{H) for J^/J^ = 3.6 and 
J' = 20 mK, the analytical mx{H) dependence defined 
by Eq. ^ is fixed. As the proportionality constant be- 
tween rrix and the measured hyperfine splitting is not 
known, in Fig. 3(c) we scale the data in vertical direction 
to match the analytical curve at the maximum near 7.5 T. 
The trend for an overall field dependence is nicely repro- 
duced except for some discrepancy between 10 T and 
13 T. However, the theoretical T = prediction is here 
compared to the data recorded at 40 mK, which may not 
be fully saturated to the zero temperature limit. Indeed, 
an extrapolation to T = [solid circles in Fig. 3(c)] from 
the temperature dependence measured above 40 mK in- 
dicates that the agreement is better than what appears 
to be. In addition, we calculate mx{H) by DMRG, treat- 
ing again the interladder coupling with J' = 20 mK in 
the mean-field approximation. As shown in Fig. 3(c), the 
result is hardly distinguishable from the analytical curve. 

The field induced 3D magnetic ordering in the gapless 
phase of BPCB below 110 mK is thus perfectly described 
in the framework of weakly coupled ID systems, treated 
as LLs. Both the phase boundary Tc{H) and the field 
variation of the order parameter mx{H) exhibit a biased 
shape. This contrasts with the dome shape of both ob- 
servables from essentially 3D spin dimer systems, such 
as TlCuCls [H] or BaCuSiaOg ^26^. Moreover, m^(i?) 
in BPCB is not proportional to Tc{H), contrary to what 
is expected for the standard mean-field description of a 
phase transition in isotropic 3D systems. This is another 
signature of the underlying ID physics. However, at low 
enough temperature, close to the critical fields the LL 
description in the gapless phase collapses as the system 
restores its full 3D character [l3| • In this region the low- 
temperature 3D magnetic ordered state is predicted to be 
a Bose-Einstein condensate 13, 27 1 leading to the linear 
niz (H) dependence [lo| , which is indeed observed in our 
experiment. 

In summary, we have investigated the experimental be- 
havior of BPCB in the gapless phase, including 3D mag- 



netic ordering occuring below 110 mK. All our experi- 
mental results are perfectly described in a single theory 
based on LLs, with no adjustment or fitting of the LL 
parameters. Magnetic field variation of T]~^, of the tran- 
sition temperature and of the order parameter (at T = 0), 
which are shown to depend substantially on the coupling 
ratio J±/J\\, provide sensitive probes for the behavior of 
the dominant spin-spin correlation function and its role 
in low temperature 3D magnetic ordering. A single set of 
LL parameters is shown to control all three observables, 
which is an essential feature of the LL model. 
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